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Reconnection-driven plasmoids in blazars: fast flares on a slow 
envelope 
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ABSTRACT 

Minute-timescale TeV flares have been observed in several blazars. The fast flaring requires 
compact regions in the jet that boost their emission towards the observer at an extreme Doppler 
factor of (5 em > 50. For TeV photons to avoid annihilation in the broad line region of PKS 
1222+216, the flares must come from large (pc) scales challenging most models proposed to 
explain them. Here we elaborate on the magnetic reconnection minijet model for the blazar 
flaring, focusing on the inherently time-dependent aspects of the process of magnetic recon- 
nection. I argue that, for the physical conditions prevailing in blazar jets, the reconnection 
layer fragments leading to the formation a large number of plasmoids. Occasionally a plas- 
moid grows to become a large, "monster" plasmoid. I show that radiation emitted from the 
reconnection event can account for the observed "envelope" of ~day-long blazar activity while 
radiation from monster plasmoids can power the fastest TeV flares. The model is applied to 
several blazars with observed fast flaring. The inferred distance of the dissipation zone is 
^diss ~ 0.3 - 1 pc. The required magnetization of the jet at this distance is modest: <x ~ a few. 

Key words: galaxies: active - BL Lacertae objects: individual: PKS 2155-304- BL Lacertae 
objects: individual: PKS 1222+216 - gamma rays: theory. 



1 INTRODUCTION 

Blazars are extremely bright and fast varying extragalactic sources 
observed throughout the electromagnetic spectrum. Depending on 
the presence or absence of observed optical emission lines they are 
classified as Flat Spectrum Radio Quasars (FSRQs) or BL Lac ob- 
jects, respectively (Urry & Padovani 1995). In either case, it is be- 
lieved that the blazar emission is powered by relativistic jets which 
emerge from supermassive black holes and beam their emission at 
our line of sight. 

Blazar variability on timescales ranging from hours to decades 
has been commonly observed at various wavelengths of the electro- 
magnetic spectrum (see, e.g., Bottcher 2007). The recent discovery 
of blazar flaring on minute timescales came as great surprise. Such 
extreme flaring has now been observed from several objects both 
BL-Lacs and FSRQs [Markarian 421 (Gaidos et al. 1996); PKS 
2155-304 (Aharonian 2007; hereafter referred to as PKS 2155), 
Markarian 501 (Albert et al. 2007; Mrk 501), and the prototype ob- 
ject of the class BL Lac (Alien et al. 2012); and the FSRQ source 
PKS 1222+216 (Aleksic et al. 2011; hereafter PKS 1222)]. Fast- 
evolving TeV flares are, therefore, a generic feature of blazar activ- 
ity. 

Ultra-fast flares pose several challenges to theoretical mod- 



els for the blazar emission. The observed minute timescale vari- 
ability is too short to originate directly from the central engine. 
Modulations of the properties of the plasma in the vicinity of the 
black hole are limited by causality arguments to be longer that the 
light crossing time of the horizon / v > R^h/c — 10 4 Mg sec, where 
W S ch = 2GM BH /c 2 and M m = 10 9 M 9 M Q . Minute-long flares are far 
more likely to originate from compact emitting regions that some- 
how form in the jet (Begelman et al. 2008; Giannios et al. 2009; 
Narayan & Piran 2012). Furthermore, for the TeV photons to es- 
cape the source in PKS 2155 the emitting blob must have a Doppler 
boosting of <5>50 towards the observer. This is much larger than the 
Lorentz factor Tj ~ 10-20 typically inferred for blazar jets from 
superluminal motions (see, e.g., Savolainen et al. 2010). Moreover, 
the TeV flaring from PKS 1222 directly constrains the location of 
the emitting region. For TeVs to escape the observed broad line 
region radiation of the source the emitting region must be located 
at scales >0.5 pc (Tavecchio et al. 2011; Nalewajko et al. 2012). 
This contraint is practically independent of the assumed geometry 
of the broad-line region (Tavecchio & Ghisellini 2012). Given the 
large dissipation distance and the large inferred energy density at 
the source, the intense flaring from PKS 1222 implies an unrealisti- 
cally large jet power unless the emitting material is, again, strongly 
boosting its emission with 6 > 50 (Nalewajko et al. 2012]]. One fi- 
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nal clue for the origin of the fast flaring is that it is observed on top 
of an envelope of longer ~day-long flares. During the fast flaring 
the flux increases by a factor of ~ a few with respect that of the 
envelope (Aharonian et al. 2007; Albert et al. 2007). 

A large number of theoretical interpretations have been put 
forward to explain the fast flares in individual sources. Fast beams 
of particles at the light cylinder (Ghisellini & Tavecchio 2008) or 
the interaction of the jet with a red giant star (Barkov et al. 2012) are 
some of them. Rarefaction waves in magnetized shells (Lyutikov 
& Lister 2010), relativistic turbulence in the jet (Narayan & Piran 
2012) or reconnection-driven minijets (Giannios et al. 2009) may 
also be responsible for the fast flares. We focus here on the latter 
possibility. 

In the MHD jet paradigm (Blandford & Payne 1982) the jet 
is expected to emerge from the central engine in the form of a 
Poynting-flux dominated flow. If the magnetic field configuration is 
appropriate after the jet acceleration and collimation phases, mag- 
netic reconnection can effectively dissipate magnetic energy in the 
flow. As pointed out in Giannios et al. (2009; 2010; Nalewajko et 
al. 2011) magnetic reconnection dissipates energy in compact re- 
gions characterized by fast motions within the jet, i.e., the radiating 
plasma can move faster than the jet material on average. The ex- 
treme Doppler boosting of the emitting region and its small size can 
naturally account for the fast-evolving flares observed in blazars. 

The reconnection minijet model is, however, based on a steady 
reconnection picture. Both observations and recent advances in 
reconnection theory reveal that reconnection is a highly time- 
dependent and dynamical process. Time dependent aspects of re- 
connection turn out to be critical in understanding the multiple ob- 
served timescales related to blazar flaring. The goal of this work it 
twofold: (i) to relax the steady-state assumptions of the reconnec- 
tion model for blazar flaring and (ii) confront the model against all 
the available observational constraints. 

In the Sec. 2 I summarize some of the recent observational 
and theoretical progress in understanding time-dependent aspects 
of magnetic reconnection. In Sec. 3 this knowledge is applied to a 
blazar jets predicting the relevant timescales and energetics of flar- 
ing. The model is applied to specific sources in Sec. 4. I conclude 
in Sec. 5. 



2 MAGNETIC RECONNECTION: A DYNAMICAL 
PROCESS 

In this Section I summarize recent progress in magnetic reconnec- 
tion theory that is relevant to the blazar jet application presented 
here. Reconnection is the process of fast release of magnetic en- 
ergy during a topological rearrangement of magnetic field lines. It 
takes place when magnetic field lines of opposite polarity are com- 
ing together towards the reconnection plane (x - y plane in fig. 1) 
and annihilate, liberating magnetic energy that heats the plasma and 
accelerates particles. The large scale / of the reconnection region is 
determined by the distance over which the magnetic field strength 
drops by a factor of ~2 (along the y direction). The magnetic pres- 
sure gradient and also magnetic tension along the y direction result 
in the bulk acceleration of the reconnected material to the Alfven 
speed Va of the upstream fluid. The fast outflow in the downstream 
allows for fresh magnetized fluid to enter the region and reconnect. 
Observationally reconnection has been extensively studied 

particles is extremely focused in the rest frame of the emitting material; see 
Nalewajko et al. (2012). 



during solar flares and in Earth's magnetotail. Laboratory exper- 
iments complement these studies in a controlled environment. A 
richness of processes that take place on very different timescales 
have been revealed by these works (see, e.g., Aschwanden 2002). A 
characteristic long timescale of the process is the global reconnec- 
tion timescale f rcc ~ I/eV A over which the magnetic energy stored 
in a region of typical scal^l / is released. Here e parametrizes the 
reconnection speed; with e ~ 0.1 been a typical observationally 
inferred value. Besides the global reconnection timescale, much 
shorter timescale variability and eruptive events are evident both 
observationally and experimentally highlighting the very dynami- 
cal nature of the process ( see e.g., Lin et al. 2005; Park et al. 2006; 
Karlicky &■ Kliem 2010) . For instance a solar flare of typical dura- 
tion of ~10 min can show strong variability on ~s timescale (see, 
e.g., Karlicky & Kliem 2010). 

For some time the reconnection theory has been dominated 
by steady-state models (Sweet 1958; Parker 1957, Petchek 1964). 
They provide intuition on how average properties such as the re- 
connection speed, the outflow speed and temperature of the recon- 
nected fluid depend on parameters. Steady state models assume 
a continuous inflow of plasma in the reconnection region and a 
smooth outflow. As such they cannot account for the erratic be- 
havior observed at the current sheet. 

A distinctly different picture has emerged from recent the- 
oretical studies of magnetic reconnection. When the resistivity 
n is sufficiently low, e.g., the corresponding Lundquist number 
S = VaIIt] » S c = 10 4 (as expected in most Astrophysical en- 
vironments) the reconnection current sheet is formally predicted 
by the Sweet-Parker theory to be extremely thin, with thickness 
d/l = S~ 1/2 <k Sc >12 ~ 0.01. Very thin current sheets suffer from 
tearing instabilities that lead to their fragmentation to a large num- 
ber of plasmoids separated by smaller current sheets (Loureiro et 
al. 2007; Lapenta 2008; Daughton et al. 2009; Loureiro et al. 2009; 
Samtaney et al. 2009; Bhattachackarjee et al. 2009; Loureiro et al. 
2012b). As a result the reconnection process is fast and resistivity- 
independent. Plasmoids grow fast through mergers and leave the 
reconnection at a speed comparable to the Alfven speed of the up- 
stream plasma (Uzdensky et al. 2010; Loureiro et al. 2012a). The 
typical plasmoid forms away from the reconnection center (the, 
so-called x-point) at y ~ I and grows up to a characteristic size 
R p <k /. Plasmoids that form fairly close to the reconnection cen- 
ter have, however, more time available to merge and grow on their 
way out of the reconnection region. These plasmoids undergo sig- 
nificant (exponential) growth to reach macroscopic scale. They are 
referred to as "monster" plasmoids (Unzdensky et al. 2010). Their 
size reaches a scale R p = S^ 1 ' 4 / - 0.1/, i.e., approaching the global 
reconnection scale. The growth timescale of the volume of a plas- 
moid is exponential e 3 ' /2,/1 , where = //Va- The mass doubling of 
the plasmoid takes place on a timescale 2 In 2t A /3 with the plasmoid 
emerging from the reconnection region at the Alfven speed Va- Re- 
sistive MHD simulations support this theoretical picture (see, e.g., 
Bhattachackarjee et al. 2009; Loureiro et al. 2012a). 

The monster plasmoids consist of energetic particles that have 
undergone acceleration in the secondary current sheets. Additional 
acceleration of particles takes place during the merging process of 
the plasmoids (Zenitani & Hoshino 2001; 2005; Drake et al. 2006, 
Sironi & Spitkovsky 2012; although the energy spectrum of the 
particles depends on the details and is still topic of investigation). 



2 For simplicity and throughout this paper I will assume that the reconnec- 
tion region has the same characteristic scale / in all directions. 
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Figure 1. Sketch of the magnetic reconnection region with a large scale /. 
Magnetic field lines of opposite polarity annihilate at the x - y plane with 
speed v rec , i.e., the reconnection proceeds on a global timescale f rec = 1/ec. 
The physical conditions relevant to AGN jets, the current layer is ex- 
pected to fragment to a large number of sub-layers separated by plasmoids 
(Loureiro et al. 2007). The plasmoids leave the reconnection region with 
the Alfven speed Va powering the "envelope" emission. Occasionally, plas- 
moids grow to become "monster" plasmoids (shaded blob) with scale ~ 0. 1 / 
giving rise to powerful, fast-evolving flares of duration ffl^e 4C ( rec . 



As I demonstrate below, the macroscopic scale of the monster plas- 
moids, their short growth timescale and fast motion and energetic 
particles that they contain make them natural candidates for pow- 
ering the ultra-fast blazar flares. 



3 APPLICATION TO BLAZAR FLARING 

In the MHD-driving paradigm for jets (e.g., Blandford & Payne 
1982), it is postulated that jets emerge from the central engine in 
the form of Poynting-fiux dominated flows. Further out the flow 
converts part of its magnetic energy into kinetic. At a (both theoret- 
ically and observationally very uncertain) distance ^diss, the blazar 
emission emerges. If the magnetic field configuration is appropri- 
ate, magnetic reconnection can effectively dissipate magnetic en- 
ergy in the flow and power the blazar emission. I assume that sub- 



stantial magnetic energy release takes place in reconnection regions 
of characteristic scale Z'.0 

The location of i? diss and the scale /' are highly model depen- 
dent. The trigger of magnetic dissipation may be instabilities that 
develop in the jet. Even if the jet is launched by an axisymmetric 
magnetic field configuration, non-axisymmetric instabilities can in- 
troduce smaller scale field structures. Current-Driven Instabilities 
(CDIs) are likely to be the most relevant in strongly magnetized 
jets (Eichler 1993; Begelman 1998; Giannios & Spruit 2006). The 
observational indications that the jet opening angle is related to the 
jet Lorentz factor through the relation OjFj ~ 0.2 (Pushkarev et al. 
2009) implies that causal contact is established in the transverse di- 
rection of a high-cr jet (where the magnetization <x = B' 2 /4npc 2 ). 
Under these conditions, CDIs can potentially grow as soon as the 
jet develops a dominant toroidal field component. CDIs are non- 
asymmetric instabilities that reorganize the field configuration. In 
the non-linear stages of their development small scale field re- 
versals may be induced in the jet allowing for energy extraction 
through reconnection (Moll 2009). The non-linear stages of CDIs 
are, however, poorly understood making it hard to predict the dis- 
tance at which they develop or the characteristic scale of the recon- 
nection layer. 

An interesting alternative is that the magnetic field is not ax- 
isymmetric at the launching region. The jet contains small scales 
field structures imprinted from the central engine (Romanova & 
Lovelace 1997, Levinson & van Putten 1997, Giannios 2010, McK- 
inney & Uzdensky 2012). Such configurations can introduce field 
reversals in the jet of the order of the size of the black hole horizon 
Rxh = 3x 10 14 Mci cm (in the lab frame). The scale of the field rever- 
sal in the rest frame of the jet is I' = r^ Sch 6 x 10 15 Tj 2 o^9 cm. 
Substantial dissipation takes place when the reconnection timescale 
free = I' Ice = ry?sch/ec becomes comparable to the expansion 
timescale of the jet f cxp = R Aiss /TjC, i.e., at a distance R iiss = 
rj^sch/f - 10 18 M 9 r^ 20 /e_i cm. I will show below that a dissi- 
pation distance of R^ ss > 10 18 cm and reconnection layers of char- 
acteristic scale of I' < 10 16 cm are within the range implied by the 
observed properties of fast blazar flares. 

3.1 Fast flares from monster plasmoids 

For the physical conditions prevailing at the reconnection region 
relevant to AGN jets, the current sheet is expected to fragment 
into a large number of plasmoids while the reconnection process 
proceeds fast (see Appendix A for details). Some plasmoids regu- 
larly grow into "monster" plasmoids, i.e., large magnetized blobs 
that contain energetic particles freshly injected by the reconnection 
process (Uzdensky et al. 2010). The relativistic motion of the plas- 
moids in the rest frame of the jet results in additional beaming of 
their emission (i.e., beyond that induced by the jet motion). When 
the layer's orientation is such that plasmoids beam their emission 
towards the observer, powerful and fast evolving flares emerge. 
Here we focus on the characteristic observed timescales and lu- 
minosities resulting from plasmoids that form in the reconnection 
region. To this end I assume that the dissipated energy is efficiently 
converted into radiation. In practice electrons are likely to be re- 
sponsible for the emission (e.g. see Nalewajko et al. 2012) so I, 
in effect, assume that a significant amount of the dissipated energy 
is deposited to electrons which undergo fast radiative cooling. The 

3 Hereafter primed quantities are measured in the rest frame of the jet while 
double-primed quantities in the rest frame of a plasmoid. 



© 2009 RAS, MNRAS 000,[TJ-?? 



4 Dimitrios Giannios 



CO 



t ~10min 



flare 




t ~1 day 

env J 



Time t 



Figure 2. Sketch of the envelope-flare structure of the emission from a re- 
connection layer. The envelope duration corresponds to that of the recon- 
nection event: f env = V /Ytec. Monster plasmoids power fast flares which 
show exponential rise and last for ffl alc = 0.1/'/<>pC. For an envelope of 
~ lday blazar flaring the model predicts that monster plasmoids result in 
~ 10-minute flares. 



latter assumption will be checked a posteriori. The former maybe 
justified by the efficient electron acceleration by the electric field 
present at the current sheets but remains an assertion. This study 
can be trivially generalized by including the efficiency factor with 
which dissipated energy converts into radiation. 

Consider a spherical blob (or plasmoid) emerging from the 
reconnection layer moving with the Alfven speed of the recon- 
nection upstream, i.e, with a corresponding bulk Lorentz factor 
y p = Vl + cr ~ a few and scale R' (both measured in the jet rest 
frame fl The growth timescale of the plasmoid t' Ax = 21n2/'/3c 
corresponds to plasmoid growth time of f^' = 21n2/'/3y p c at its 
own rest frame. The observed rise time for the plasmoid emis- 
sion is frise.obs — KLc^p ~ 0.51' /y p 6 p c, where S p is the Doppler 
boost of the plasmoid radiation towards the observer. The plasmoid 
subsequently cools down and expands on an observed timescale 
'deciine.obs ~ R p l c 5 p - For the typical the typical parameters y p ~ 2, 
R' p = 0.1/' adopted in what follows, the observed rise and decline 



4 I treat the plasmoid as a sphere in the jet rest frame instead of it own rest 
frame. It is unclear which approximation describes better the reconnection 
plasmoids for relativistic reconnection. In the limit of modest relativistic 
motions (y p ~ a few), this distinction is not affecting the results presented 
here by much. Hereafter I assume that R' p 



timescales are similar. For simplicity, we define as the variability 
timescale t v = Wimcobs = R%/c6 p . 

Consider a jet emerging from a supermassive black hole with 
luminosity Lj, opening angle dj and Lorentz factor r,-, We also as- 
sume that OjTj = 0.2 as indicated by observations (Pushkarev et al. 
2009). The typical bulk Lorentz factor of gamma-ray active blazars 
is Tj ~ 10 - 20 (Savolainen et al. 2010; Piner et al. 2012). The 
energy density at the dissipation, or "blazar", zone is 



U' = 

J ne 2 r 2 R^c 

J J 



: 0.26^ 



erg/s, 



(1) 



where the Lj is normalized to 10 46 erg/s and the distance R to pc in 
anticipation of the limits set by observations. 

Pressure balance across the reconnection layer requires the 
energy density of the plasmoid to be similar to that of the jet 
Up ~ U'^ Assuming efficient conversion of dissipated energy 
into radiation, the rest-frame luminosity of the plasmoid is thus 
L" = U' p '(4n/3)R p 2 c. This luminosity can be converted to the 
observed luminosity Lp, hs = 8 p L". Because of the R' p 2 depen- 
dence of the luminosity it is clear that the largest "monster" plas- 
moids (with R'^ = 0.1/') power the brightest flares. The size of 
the plasmoids is directly constrained from the observed variabil- 
ity timescale R' p = t v S p c. Putting everything together, the observed 
luminosity of the plasmoid is 



->,50'i 



-2 
'PC 



erg/s, 



(2) 



where 6 p = 50(?5o and R = li? pc pc, and t v = 10f,,_io min. 

The Doppler factor 6 p depends on several parameters. It is re- 
lated to Tj, y p , the angle of the plasmoid with respect to the jet 
motion and the observer's angle of sight. For perfectly aligned jet, 
plasmoid and observer S p ^ 4rjy p . In, perhaps, more common 
situations where the reconnection layer is at a large 9 ~ n/2 an- 
gle with respect to the jet propagation (as seen in the jet rest) and 
fairly aligned with the observer (giving powerful flares) 6 p ~ Tjy v . 
For the demonstrative arithmetic examples used here we adopt 
S p = 1.25rjy p = 50rj_2o7 P .2- One can see (see eq. 2) that powerful 
flares on a timescale of ~10 min is possible even with very modest 
relativistic motions within the je0 y p ~ 2. 



3.1.1 Ejection of multiple monster plasmoids 

During a reconnection event multiple monster plasmoids are ex- 
pected to form. The seed of a monster plasmoid forms fairly close 
to the reconnection center region y <k /' and spends sufficient time 
in the reconnection region to grow to a large size. 2D simulations 
(Loureiro et al. 2012a) indicate that monster plasmoids form ev- 
ery a few Alfven times t A or at a rate of ~ 0.3^'. It appears likely 
that 2D simulations underestimate the rate of formation of monster 
plasmoids. The actual rate may be higher when the 3D structure of 



5 The exact relation of Up and U'. depends on the magnetic and heat con- 
tent of the upstream and the plasmoid. For instance balancing the pressure 
of cold, strongly magnetized upstream pj = B' 2 /87T = U'./2 with that of 
an assumed relativistically hot plasmoid p p = U p /4, we find U' p = 211'-. 
The expression is slightly different if the magnetic field contributes to the 
pressure of the plasmoid. The Assumption U' p ~ V, is expected to hold 
within a factor of ~ 2 independently of these details. 

6 Narayan & Piran (2012) come to a different conclusion concerning the 
fastest flares expected by magnetic reconnection. Their analysis assumes 
steady reconnection and is therefore applicable to the global reconnection 
timescale (topic of next Section) but is not relevant for the fast flares pow- 
ered by plasmoids. 
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the layer is considered. The x-point is in reality an elongated struc- 
ture (along the x-axis of fig. 1) providing a larger physical region 
where the seeds of monster plasmoids form. Monster plasmoids po- 
tentially form at a rate up to I' /R' ~ 10 higher than that found in 
2D studies. Clearly this question can only be answered high resolu- 
tion, 3D restive MHD simulations. If monster plasmoids emerge at 
arate ~ (0.3-3);^', some (3-30)/e_i plasmoids are expected from 
a single reconnection layer powering multiple flares. The observed 
properties of the monster plasmoids are determined by the basic 
properties of the reconnection region that generates them. To the 
extent that all monster plasmoids reach similar size of ~ 0.1/', the 
model predicts a similar duration and brightness for this sequence 
of fast flares. A sketch of such pattern is shown in Fig. (2). 

3.2 The "envelope emission" from the reconnection region 

The bulk motion of a monster plasmoid is expected to be similar to 
the speed of other structures (e.g. smaller plasmoids) leaving the re- 
connection region. When the plasmoid emission is beamed towards 
the observer (powering a fast flare), the overall emission from the 
current layer is also beamed by a similar factor. The emission from 
the layer forms a slower-evolving envelope. In the following we 
calculate the timescale and luminosity of the emission from the re- 
connection layer. 

The size of the reconnection region is constrained to be 
/' = 10/?p =a 10 16 f v ,iorj,2o7 P ,2 cm. The observed time it takes for 
the reconnection process to be completed is t Ixobs = /' /Tjec = 
1.5 x 10 5 y p 2f v ,io/ f -i s - F°r fast flares lasting t v ~ minutes, the en- 
velope emission is naturally expected to last hours or days. If the 
reconnection timescale is longer than the expansion time of the jet 
fcxp.obs — ^diss/2r^c, the reconnection is likely to be interrupted by 
the global expansion of the plasma. In this case ? C xp,obs determines 
the observed duration of the eventQ The competition of the recon- 
nection and expansion timescales of the jet set the duration of the 
envelope: 

t em =min - — , — r- = mm 1.5X10 3 — , l.SxlO 5 -;— s.(3) 

lijec TTjci I El T z j2Q \ 

The duration of the envelope emission lasts for ~days. Such 
timescale is characteristic of blazar flares (See next Section). 

The (lab frame) energy available to power the envelope emis- 
sion is £ c „v = U s 2l' 3 IT s , where Uj = T 2 ,U'j is the energy density 
of the jet and 2/' 3 /Ti accounts for (lab frame) volume of the recon- 
nection region that powers each minijet (see fig. 1). The emitted 
luminosity of the reconnection region is E tm /t ewl . It can be con- 
verted into observed luminosity by accounting for beaming factor 
of the emission 6 2 p : 

L cm , obs = lt\6\l a V\ec = LSxlO^r^Xo^iof-i^^o erg/s.(4) 

The envelope emission is quite bright. Dividing eqs. {2]l and 
l(4j, one arrives to a very simple expression for the ratio of the plas- 
moid to envelope luminosities L p /L em ^ ly^/f-i- The luminosity 
contrast depends only on the Lorentz factor of the minijet in the 
rest frame of the jet and the reconnection rate e. As we discuss in 
the next Sections, this ratio is observed to be of order unity con- 
straining 7 P ~ a few for e ~ 0.1. 



So far I have considered the emission from a single recon- 
nection region which is beaming its emission towards the observer. 
When the reconnection scale is smaller than the cross section of the 
jet (/' << R6j), there may be as many as ~ (Rdj/l') 3 reconnection 
regions in the emitting zone. Even after corrected for beaming of 
the emission from current sheets, up to ~ (R6j/l') 3 /4yj; reconnec- 
tion regions may beam their emission towards the observer. The 
overall amplitude of the envelope emission and the number of fast 
flares are therefore enhanced by this factofj In this case, the con- 
trast L p ,obs/ienv,obs drops since more than one reconnection layers 
contribute to the envelope emission. 



4 APPLICATION TO OBSERVATIONS 

In this section we examine how the wealth of information of fast 
TeV flaring among blazars can be used to extract information on 
the physical conditions of the emitting region and constrain the re- 
connection model. We first discuss observations of several BL Lac 
objects and then FSRQ PKS 1222. 



4.1 Flaring BL Lacs 

Fast TeV flares have been observed in several BL Lac objects in- 
cluding PKS 2155, Mrk 501, Mrk 421 and BL Lac itself. The vari- 
ability timescale ranges from a few to ~10 minutes. Most detailed 
are the observations of PKS 2155 which reveal an envelope emis- 
sion of duration of hours that contains several ~ 5-minute flares 
of comparable luminosity (Aharonian et al. 2007). Mrk 501 also 
shows flaring with similar flare-to-envelope ratio. BL Lac, the pro- 
totype source has also been observed to flare on 10 minutes sug- 
gesting that minute-flaring is a generic property of blazars. 

First we apply the model to the PKS 2155 observations (see 
Aharonian et al. 2007). The observed luminosity of the envelope 
and of the fast flares is Lj. nv ~ L ff ~ 3x 10 46 erg/s. The fast flares last 
for t v ~ 5 min. When the observation started the envelope emission 
was already on and lasted for ~100 min. Narayan & Piran (2012) 
use Bayesian analysis to estimate the mean duration of the enve- 
lope emission of f env ~ 2 x 10 4 sec. The true jet power is uncertain 
but Lj ~ 10 45 erg/s appears a reasonable estimate given observed 
luminosity of the source of up to 10 47 erg/s (to be corrected for rel- 
ativistic beaming). A high Doppler factor 6 V > 50 of the emitting 
material is required for the escape of the TeV radiation from the 
soft radiation field of the jet without extensive pair creation (Begel- 
man et al. 2008). Observations do not directly constrain where the 
emission takes place in this source. 

The fact that L cnv ~ L s means that y p ~ 2 (compare eqs. (2) 
and (4)). Setting S p = 50 and the observed f v and inferred L s in eq. 
(2), we derive R^ ~ 7x 10 17 cm. With the parameters of the model 
determined, the envelope duration is predicted to be f cnv ~ 3 x 10 4 
sec (see eq. 3). This timescale is very close to the one inferred from 
observations. PKS 2155 showed (i) several fast flaring events of (ii) 
similar characteristic timescale and luminosity. Multiple flares have 
a natural explanation within the reconnection model. They can be 
understood to come from different monster plasmoids that emerge 
from the same reconnection region. 



7 The reconnection process takes about an Alfven crossing time to establish 8 The opposite limit where /' > R8j is not physical since the reconnection 

itself. For yet faster expansion ? exp < = frecM the reconnection process region should fit within the jet cross section: the condition I' < R9j must be 

may not have sufficient time to develop in the first place. satisfied. 
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4.2 TeV flares from PKS 1222 

The model can also be applied to the FSRQ PKS 1222. In this 
source the dissipation distance is robustly constrained to be Ra ss > 
0.5 pc (Tavecchio et al. 2011; Nalewajko et al. 2012). If the emit- 
ting region is characterized by 5 ~ Tj ~ 20 the required luminosity 
of the jet is unrealistically high (Nalewajko et al. 2012). A moder- 
ately higher S ~ 50 (in agreement to that inferred for PKS 2155) 
is, however, sufficient to relax the energetic requirements of the jet 
and is adopted here. During the epoch of the TeV flare there is 
an envelope of high y-ray activity. Fermi-LAT detected a flare of 
L cnv ~ 10 48 erg/s and duration of roughly ? cnv ~ 10 5 sec (Tanaka 
et a. 2011). The fast flares are observed with MAGIC in the sub- 
TeV energy range: Lg ~ 10 47 erg/s (with total luminosity possibly 
higher by ~ a few to account for bolometric corrections in view of 
the steep observed TeV spectrum). The flux evolved on a timescale 
of t v ~ 7 min. The true jet power is also uncertain but Lj ~ 10 4S 
erg/s appears a conservative estimate given that the observed disk 
emission been a factor of several larger and that the beamed lumi- 
nosity of the source reaches 10 48 erg/s. 

Setting (5 p = 50, L s ~ 10 46 erg/s and f v ~ 7 min in eq. (2), 
the dissipation distance is found to be R iiss ~ 1 pc. The envelope 
emission is predicted to last for (see eq. (3)) ? env ~ 10 5 sec that 
is consistent with the observations. The observed envelope emis- 
sion is observed to be more luminous than the fast flare by a factor 
of several (though weaker, the fast flare is clearly observed with 
MAGIC because of its harder emission). Given the adopted param- 
eters, eq. (4) implies that a single reconnection region has envelope 
luminosity L cnv ~ Lff. Possibly several reconnection layers con- 
tribute to the envelope emission simultaneously since ROj/l' ~ 5 
enhancing the ratio of the luminosity of the envelope emission with 
respect to that of fast flares. 

Summarizing all blazar flares can be accounted for by little 
changes of the physical parameters of the model. Typically I infer 
^diss ~ (0.3 - 1) pc, 7 P ~ 2 and the size of the reconnection region 
Z'<10 16 cm. In the Sect. 5 we discuss, among other issues, a possible 
physical motivation for these values of the parameters. 

4.3 Radiative mechanisms and particle cooling 

The energetic requirements for the fast flaring can become more 
stringent if the radiating particles (assumed to be electrons) are 
not in the fast cooling regime. Here we assume that the TeV emis- 
sion is either result of synchrotron self Compton (SSC) or external 
inverse Compton (EIC) and investigate the electron energetics re- 
quired to produce the observed ~ lOOGeV- multi TeV emission and 
whether they are likely to radiate efficiently for the model parame- 
ters adopted in the previous Section. 

To assess whether the emitting particles "cool fast", the ex- 
pansion timescale of the plasmoid t" xp = Rp/c = 6 p t v = 3 x 
10 4 cS5of v ,io s is to be compared with the radiative cooling timescale. 
In the case of SSC emission from e~ (or pairs) with random 
Lorentz factor y e in magnetic field B", the characteristic en- 
ergy is esse - 10~ 8 (5y 4 B" eV. Depending on details of the re- 
connection configuration (such as guide field strength), the plas- 
moid can be magnetically or heat dominated^ For simplicity, I 
parametrize the magnetic field strength in the plasmoid as B" = 

9 The fact that the jet is Poynting-flux dominated in this model does not 
mean that the emitting region is also strongly magnetized. On the contrary, 
efficient reconnection may result in weakly-magnetized (heat-dominated) 
plasmoids. 



(e B 47rt/;') 1/2 = l^ /3 l}£R£ Gauss. Setting e ssc = 100e 100 GeV 
GeV, one finds for the required electron (random) Lorentz factor 
y e = 2 x lO 4 ^!^/?^ 4 e^ 8 L 46 1/8 5j 1/4 . The SSC cooling time scale 
is q &c = 5 x 10 8 /(1 + y)y e B' 2 s or 

3 3 Rl'Xo 
fs'sc - § x 10 — -775 — -7/8—7/8 s ' © 

J e 100GcV fc B,l/3 i -'46 

where y ~ a few accounts for the ratio of the SSC to synchrotron 
power. 

If a substantial external radiation field of energy density C/ rad 
is present, it can contribute to the particle cooling through EIC. 
Assuming an isotropic radiation field of characteristic energy e see d, 
the energy of the up-scattered photon is e Elc — <5pTe e seed (f° r scat- 
tering in the Thomson limit). Solving for the electron Lorentz fac- 
tor: y e = 7 x 10 3 (emc.iooGcv/<? S!: cd,i) 1/2 <55o, where e secd = le seed ,i 
eV. The energy density of radiation in the rest frame of the blob 
ls k'rad ~ S 2 Ut*a- The EIC cooling timescale for such electron is 
'mc = 2xl0 7 / rf (/ r '; d sor 

i.4xio 4 / eseKU y/ 2 

hie = 77 T\- s > (6) 

<Aeed,-4O50 V ^lOOGeV ' 
for (/ S ecd = 10~ 4 t/ sced ,-4 erg/cm 3 . 

In the case of PKS 2155 no powerful ambient isotropic ra- 
diation field is evident. The plasmoid may well propagate into a 
dense radiation field emerging from the large scale jet or other re- 
connection regions. This depends however on uncertain details of 
the overall geometry (see Nalewajko et al. 2011 for various possi- 
ble geometrical configurations). On the other hand SSC emission 
has to be present. Setting i? diss ~ 7 x 10 17 cm, Lj ~ 10 45 erg/s (see 
previous Section), and e B = 1/3, y = 2, I arrive at t ssc ~ 5 X 10 3 
s (see eq. 5). This timescale is, by a modest factor, shorter than the 
expansion time of the blob f" xp = 1.5 x 10 4 £50^,5 s indicating that 
efficient ~TeV emission is plausible. The required random Lorentz 
factor of the emitting electrons is y e > 10 4 . 

One can derive a similar estimate for the SSC cooling 
timescale for the parameters relevant to PKS 1222 flaring, i.e. 
f/ssc < f" xp ~ 2 x 10 4 s. Another obvious possibility for emission 
from this source is EIC of photons from the infrared torus. With 
U/r ~ 10~ 4 erg/s and characteristic energy ei R ~ 0.3 eV, the EIC 
takes place in the Thomson regime (Nalewajko et al. 2012). The 
typical electron emitting ~ 100 GeV has y e ~ 10 4 with a cooling 
time scale of the particles of t'^ ic ~ 8 x 10 3 s ~ t'^ sc < f" xp . 

In Summary the ~TeV emitting electrons are characterized 
by a random y e ~ 10 4 and have a cooling timescale somewhat 
shorter than the expansion timescale of the blob allowing for ef- 
ficient TeV emission. The detailed spectrum necessarily depends 
on the details of the particle distribution that are model depen- 
dent. However, an equipartition argument (e.g. sharing of the dissi- 
pated magnetic energy between electrons and protons) would give 
y e ~ (m p /2m e )<T = 3 x 10 3 o" 3 for the electrons in the plasmoid, 
where a = 3<r 3 is the upstream magnetization. Therefore, a modest 
particle acceleration above equipartition is sufficient to explain the 
observed emission. 



5 DISCUSSION/CONCLUSIONS 

The similarities in the variability patterns seen in several blazars 
(PKS 1221, PKS 2155, BL Lac, Mrk 501) are striking: fast TeV 
flares on -minutes timescale that appear on top of an envelope of 
enhanced gamma-ray activity that lasts for hours or days. The simi- 
larities strongly indicate similar physical conditions at the emitting 
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region: large Doppler factor 5 P ~ 50 and a dissipation zone located 
at ~pc distance from the black hole. 



5.1 Models for fast blazar flaring 

Several suggestions have been put forward for the ultrafast blazar 
flaring. Fast electron beams (with bulk y e ~ 10 6 ) may develop 
along magnetic field lines close to the light cylinder (i.e. within sev- 
eral Rsa, see , Tavecchio & Ghisellini 2008). The TeV flare, in this 
scenario, is result of the beam inverse Compton scattering external 
radiation fields. This model fails to account for the large emission 
distance required by PKS 1222 (Aleksic et al. 201 1). Alternatively, 
a red giant may cross the jet (Barkov et al. 2012). The ram pres- 
sure of the jet strips the envelope of the star which consequently 
fragments. Emission from shocked stellar and/or jet plasma may 
power blazar flares. While stellar-jet encounters are expected, any 
interaction region will necessarily move slower that the jet r inl <T, . 
The required Doppler boost of the emitting region towards the ob- 
server S ~ 50 may therefore be hard to explain in this picture. Al- 
ternatively, a recollimation shock on pc scales can help to focus 
the jet energy in a small region inducing short timescale viability. 
However, non-axisymmetries in the jet-external medium interac- 
tion will likely make the focusing insufficient to explain the most 
extreme flares (Bromberg & Levinson 2009). If the jet activity is 
sufficiently erratic, the jet can be envisioned as magnetized shells 
separated by vacuum (Luytikov & Lister 2010). Rarefaction waves 
of the front part of a shell can reach a bulk Lorentz factor much 
higher than that of the jet on average. Fast flares may come from 
these rarefied regions. Relativistic turbulence in the jet can also al- 
low for emitters-blobs moving with T b » Tj to be responsible 
for intense and fast flares (Narayan & Piran 2012). For the turbu- 
lence not to be supersonic (or it would decay fast by shocks) the jet 
must be Poynting flux dominated. The driver and the region where 
the turbulence develops remain to be identified. Magnetic recon- 
nection could drive the turbulent motions (and turbulence can, in 
turn, enhance the reconnection rate; Lazarian & Vishniac 1999). In 
this case, however, it is quite likely that the most powerful flares 
are directly related to the driver, i.e., the reconnection process itself 
rather than the turbulent motions. 



5.2 The reconnection model for fast flaring 

Here we have revisited the reconnection minijet model for the fast 
flaring (Giannios et al. 2009; 2010). We focus on time-dependent 
aspects that are naturally expected (and directly observed in labo- 
ratory experiments and solar system environment) in the reconnec- 
tion process. It is demonstrated that at least two timescales appear 
in the problem. The longer one is associated with the time it takes 
for a magnetic energy to be dissipated in the reconnection region 
and creates an "envelope" of flaring activity that lasts for several 
hours to days. Instabilities in the current sheet (e.g. secondary tear- 
ing instability; Loureiro et al. 2007) result in erratic formation of 
plasmoids that leave the reconnection region at relativistic speed. 
The largest ones, "monster" plasmoids, can power the fast, minute- 
long blazar flares. Several to tens of monster plasmoids can emerge 
from a single reconnection layer. The super fast flaring may there- 
fore not happen in isolation. A sequence of fast flares are expected 
to have similar timescale set by the size of the reconnection layer 
as observed in PKS 2155. 

The model can be applied to all blazar sources for similar 
adopted parameters. Dissipation can take place at distance R& ss ~ 



0.3 - 1 pc, the bulk Lorentz factor of the jet is Tj ~ 20, and the 
size of the reconnection region /' < 10 16 cm. These quantities point 
towards an interesting possibility for the magnetic field structure 
in the jet and the origin of the blazar emission. If the magnetic 
field configuration is not exactly axisymmetric at the horizon, the 
jet may emerge with small-scale field structures of size similar to 
that of the central engine ~ i?sdi ~ 3 x 10 14 M9 cm (along the direc- 
tion of the jet propagation). Even a modest non-axisymmetry at the 
base of the jet can be amplified by stretching in the toroidal direc- 
tion because of the jet (lateral) expansion. The jet expands from a 
lateral size of r ~ Rs,:h at the launching radius to r = 8jR^ ss » ^sch 
at the dissipation distance. The resulting scale of the field reversals 
in the rest frame of the jet is ~ T^sch that may be used as an esti- 
mate of the characteristic scale of the reconnection layer /'< 10 16 cm. 
For typical parameters, the reconnection time catches up with the 
expansion time of the jet at distance R^ ss ~ T^Rs^Je ~ lpc. 

The bulk motion of plasmoids in the jet rest frame required 
for the model to work are very modest y p ~ 2. The bulk mo- 
tion in the reconnection picture corresponds to the Alfven speed: 
7 P - Vl + cr, implying that a magnetization cr ~ 3 is required at 
the dissipation zone R^ B . Is it reasonable that the jet remains mod- 
estly Poynting-flux dominated at ~ lpc scale? This would imply 
that the conversion of Poynting-to-kinetic flux is not complete by 
that distance. A systematic study of superluminal motions on pc to 
tens of pc scales reveals that blazar jets still undergo acceleration 
out to the larger scale (Piner et al. 2012). In the context of MHD 
acceleration of jets, this would imply that indeed the pc-scale jet 
maintains a substantial magnetization. 

Ultra-fast flares are the tip of the iceberg of blazar variabil- 
ity. The process of magnetic reconnection is potentially responsi- 
ble for powering a much broader range of blazar activity. Recon- 
nection may well take place at larger (e.g. multi pc) scale where 
the plasma is presumably less magnetized (because of further con- 
version of magnetic energy into kinetic). When the reconnecting 
plasma is characterized by cr < 1, the reconnection speed v rec slows 
down (since v rec oc V A < c). In this case, the reconnection timescale 
becomes longer and reconnection layers may power days-to-weeks 
long flares of "envelope" emission. 

5.2.1 Other Astrophysical implications 

These ideas of plasmoid-dominated reconnection may be applied to 
other Astrophysical sources. The variability patterns of gamma-ray 
burst (GRB) emission show fast flaring on top of slower-evolving 
lightcurves that may be connected to such reconnection process 
(see also Lyutikov 2006; Lazar et al. 2009; McKinney & Uzden- 
sky 2012; Zhang & Yan 2012). Similar considerations may apply 
to flares observed during the GRB afterglows (Giannios 2006). Re- 
connection minijets may also be the key to understand the fast GeV 
flaring of the pulsar wind nebula of Crab (Clausen-Brown & Lyu- 
tikov 2012; Cerutti et al. 2012). In particular such model can at- 
tempt to explain the day-long flares and shorter timescale variabil- 
ity observed during major flaring of the Crab (Abdo et al. 2011; 
Tavani et al. 2011; Buehler et al. 2012). 

5.2.2 Open issues 

This study focused on the rough energetics and timescales of 
plasmoid-dominated reconnection in blazar jets. While the feasi- 
bility of the process to account for blazar flares has been made, a 
more detailed comparison to observations requires progress in our 
theoretical understanding on a number of fronts. 
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Where is the dissipation zone of jets located? Studies of the 
global jet structure and stability can reveal where and why recon- 
nection in the jet develops. These studies will also probe the char- 
acteristic length scales and orientation of the reconnection regions. 
The plasmoid-dominated reconnection is also a study in progress. 
Better understanding of fragmentation instabilities of the current 
sheet requires high-resolution 3D simulations. The theory should 
be tested and extended to the, interesting here, trans-relativistic 
<r ~a few regime. Finally for making predictions on the spectra of 
the resulting emission and direct comparison to observations, a bet- 
ter understanding of particle acceleration in reconnection regions 
is required. Particle-in-cell simulations make rapid progress in this 
direction (Zenitani & Hoshino 2005; Drake et al. 2006; Sironi & 
Spitkovsky 2011). 
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APPENDIX A: PROPERTIES OF THE RECONNECTION 
REGION IN BLAZARS 

Here I estimate the Lundquist number 5 = I'vaIt] that characterizes 
the reconnection region in a blazar jet as well as several physical 
scales of relevance to the structure of the current sheet. I focus on 
modestly high-cr upstream (i.e., a ~ a few). The resulting Alfven 
speed v A = c. The size of the reconnection region is constrained 
from the observed duration of the blazar flares to /' < 10 16 cm. 

The resistivity rj contains contributions from Coulomb colli- 
sions r} s (Spitzer resistivity) and electron scattering by photons rjc 
(the, so called, "Compton drag"; Goodman & Uzdensky 2008). For 
any reasonable temperature of the electrons in the jet, the Compton 
drag is expected to dominate over the Spitzer resistivity. Follow- 
ing Goodman & Uzdensky (2008) rjc = ccTU' mll /3ne 2 n c . For the 
very rough, order-of-magnitude estimate performed here we can 
set U' mi ~ U' (i.e., much of the energy density of the jet -as given 
in eq. (1)- converts to radiation) and n e = n p = U'-/(l + cr)m p c 2 
for an electron-proton jet. For R iiis ~ lpc the resulting resistivity 
rjc ~ 10 5 cm 2 /s. The Lundquist number is S ~ 10 21 clearly far 
in excess the critical value S c - 10 4 above which the secondary 
tearing instability in the current sheet sets in (and S » 10 5 as re- 
quired for monster plasmoids to develop). Reconnection is blazars 
is, therefore, likely to take place in the plasmoid-dominated regime. 

As discussed in Uzdensky et al. (2010) the plasmoids have 
a hierarchical structure with a large range in size. Large plas- 
moids are separated by a secondary reconnection layer of length 
Z <2) where, in turn, smaller plasmoids form. This hierarchy repeats 
on ever smaller scales / (,,) . The smallest structure would formally 
be a "critical layer" of length l c = S c r]/c ~ 1 mm with nominal 
thickness S c = S^ l! 'l c ~ 10~ 3 cm! Such small structures are never 
realized in practice because the characteristic plasma scales (e.g. 
the ion skin depth dj) are much larger. For the expected plasma 
density of the jet dj ~ 10 s cm. Once the thickness of a sublayer 
5 ~ d,, the resistive MHD description fails. One instead deals 
with "collisionless" reconnection where the resistivity is likely to 
be controlled by plasma instabilities. In this regime, the Petschek 
model for (relativistic) reconnection (Lyubarsky 2005) may ap- 
ply (Kulsrud 2001). The Petschek reconnection rate at a layer of 
length / = S' l2 dj ~ 10 10 cm is v rcc ~ c/ln5 p ~ 0.03c, where 
S p = Ic/r] ~ 10 15 . Moreover, when 8 ~ dj, Hall MHD terms be- 
come important further increasing the reconnection rate (see, e.g., 
Malyshkin 2008). This local fast reconnection rate controls also the 
overall rate of large-scale reconnection (Uzdensky et al. 2010). 
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